#ifndef __Sound_Wave_HH__
#define __Sound_Wave_HH__
#include "../../includes/equation-base.hh"

class SoundWave : public EquationBase {
public:
  SoundWave();
  virtual ~SoundWave();  
  void get_equation_flux(int dim, double* dir, double* var, double* fequ);
  // void get_local_vmax(int dim, double* dir, double* var, double& a_local_vmax);
  void get_local_vmax(int index);    
private:
  void get_1d_fx(double* fequ, double* var);
  void get_2d_fx(double* fequ, double* var);
  void get_2d_fy(double* fequ, double* var);
  void get_3d_fx(double* fequ, double* var);
  void get_3d_fy(double* fequ, double* var);
  void get_3d_fz(double* fequ, double* var);        
};
#endif // __Sound_Wave_HH__
